Green's Function of Anyons in Calogero Model and Quantum Hydrodynamics 
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We find that correlation functions at one dimension are crucially affected by the curvature of 
the bare single particle spectrum. Parabolic curvature leads to two closely related phenomena: 
the Green's function exhibits oscillation (as a function of the coordinate), while the polarization 
operator acquires support in part of the frequency-momentum plane. We calculated the Green's 
function using the WKB approximation for collective variables theory Q. Within this approach, 
the single particle Green's function is related to a quantum soliton 0]. The finite support of the 
polarization operator is due to periodic density waves. 



Tomonaga-Luttinger liquid has been a paradigm of non 
Fermi liquid systems. The bosonization technique has 
been a particularly efficient way to study it. Recently it 
was pointed out [2j that for some one dimensional phe- 
nomena the finite curvature of conductance band is es- 
sential. Such are photo-voltaic effect, thermo-power low 
temperature Coulomb drag etc. The violation of particle- 
hole symmetry is necessary for this and alike effect to 
exist. Since conventional bosonization relies on the lin- 
earization of the single-particle spectrum such phenom- 
ena can hardly be approached. 

In this work we address the issue of spectrum cur- 
vature for correlation functions of Calogero-Sutherland 
(CS) model 4]. This model appears in various branches of 
hysics, such as spin chains, FQHE, disordered metalsj^, 
The Hamiltonian of CS is 
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The special values of it A = 1,2,4 corresponds to the 
Dyson's symmetry classes of Random Matrix Theory 
(RMT) 7]. The exact eigenstates and eigenfunctions 
of CS model have been found, first by Forrester^] for 
integer and later on by HaQ for rational values of cou- 
pling strength A. The continuous description of correla- 
tion function for generic values of A is absent. A certain 
progress had been reached in the Ref.JTjj based on the 
replica symmetry breaking technique. 

In this work we focus only on long range properties 
of CS model. The structure of this paper is following. 
Following Awata 01 we P ass to the variables 1]. We 
explain what part of correlation functions we are after 
on the example of free fcrmions. We calculate Green's 
function for the large value of the interaction strength 
using saddle-point approximation. Finally we determine 
the support of the polarization operator based on the 
hydrodynamical calculations. 

We consider a ring geometry of the system. The co- 
ordinates of the particles can be described by a complex 
variable z n = Le " where 9 is an angle along the circle 



of radius L. The Hamiltonian can be rewritten as 
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The excited states = f o^f , where $^ is a symmetric 
function of the particle coordinates parameterized by 
partition k. The Hamiltonian H = 5 , ( ^" 1 _ff\E'o is acting in 
the space of symmetric (bosonic) function. After some 
algebra one arrives to 
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where D t — z^di- To prepare for second quan- 
tization one defines so called collective variables 
V{0) = Eti 5 (° ~ 0i\ Pk = JT d9e ik9 p(9). In terms of 
collective variables the Hamiltonian has a form 

d 2 



mnp n+m 



1 N 

h=- y 

2 ^ 

m.n— — N 
^ AT— 1 N—m 

m— n— 1 



Opndpn 



N 

+(i-a) y; n 2 Pn 



dp n 



PnPr, 



Opn-I 



+ p^ n p- m 



dp-n- 



(4) 



So far our calculation have been exact. Now we pass to 
the hydrodynamic limit. In the limit of N — > oo the 
Hamiltonian eq.|@J can be written as 



H = dx 
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Here x is a coordinate along the circle (x = 

p(6) etc. The 
potential energy of quantum liquid is given by 



j-9), the linear density p(a 
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transform is defined as p H (x) = —P f dx' p ^ x ; . The 
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modes of velocity operator are defined as v„ 
2tt ^— n a ^_ - + |p n sgn(n)Y It is easy to see that with 
these definitions particle velocity and density satisfy a 
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standard commutation relations one expect for any hy- 
drodynamic theory ^3] [ v i x )i p(y)] — ~~ iS'(x — y). To 
reproduce the conventional bosonization theory one ex- 
pands the Hamiltonian up to second order in fluctuations 
(due to conservation laws the velocity v and density op- 
erators Sp are of the same order) and spatial gradient 
(keeping only the terms that have no gradient). By do- 
ing so we see that on the conventional Luttinger model 
interaction enters only through sound velocity of bosonic 
modes. 

Now we use the Hamiltonian J5J) to study effects of 
the spectrum curvature. We start with single particle 
Green's function G(x, t) = -i(T^f A (x, t)^ A (0, 0)), where 
"J is an anyonic (for A =/= 1) field. Interesting on its own 
(it's spectrum function can be measured in angle resolved 
tunneling experiments (l3jl Green's function will also be 
helpful for study of more involved correlators. As we will 
see shortly, the Green's function for CS model and for free 
fcrmions turn out somewhat similar. We consider non 
interacting case first. By the Wick theorem the Green's 
function can be easily evaluated For parabolic dispersion 
the time ordered Green's function (t > 0) is given by 

G = -^^e"/ 4 e ^+^(^*) 2 Erfc ^^^(x+vt) 2 ) 

The Green's function is a sum of left and right moving 
terms, multiplied by fast oscillating term. The smooth 
part of the Green's function for right moving particles 
is plotted on Fig(l). The left tail of Green's function is 
captured by linear bosonization. There is no right tail, 
because we deal with strictly parabolic spectrum. The 
oscillations are universal. The right tail will appear as 
soon as negative spectrum curvature is take into account, 
to match a linear bosonization. A transition between 
oscillating a power law region is model dependent. 

One can get a good idea about function using a saddle 
point approximation for the integral over momentum. 



G=-; 



e lpFX e l{x ^" t)2 /2t 6{x-vt) + (L «-> R) 



One sees that saddle point gives a dominant contribution 
only if the particle moves faster than the Fermi velocity. 
In this region the result is oscillatory function. For the 
particles moving slower than Fermi velocity the saddle 
point lies out of the region of integration and the inte- 
gral is determined by its low limit. In that case we restore 
a standard power decay of the Green's function. The re- 
sult for particles moving vaster than Fermi velocity par- 
ticle does not agree with linear bosonization prediction 
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The reason of this 



disagreement is simple. We have evaluated the Green's 
function of free fermion with a parabolic spectrum. In 



metals electrons are subject to periodic potential. They 
form Bloch states and have periodic spectrum that has 
a negative curvature close to the boundary of Brillouin 
zone. This limits the maximally accessible velocity of the 
classical trajectory. The standard bosonization considers 
only the situation in which there is no classical trajectory 
available. 

Indeed, within a tight-binding model the single parti- 
cle spectrum is given by e(p) = W (l — cos (§f )) > where 
a is a lattice constant and W is a band width. We see 
that integral has a stationary point at - 
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The right hand side of this equation is a periodic function 
of momentum and it is smaller than one. So the saddle 
point contribution disappears if the velocity on the classi- 
cal trajectory exceeds the one available by band structure 
(haW -C (x — vpt)/t). This is a linear bosonization limit, 
intensively studied in the literature so far 003. We 
will consider the situation when the classical trajectory 
dominates the result, but the distance to the light cone 
is still large (vp <C (x — Vpt)/t <C fiaW). In this case 
the semiclassical approximation is justified. Close to the 
light cone ((x — vpt)/t <C vf) WKB approximation fails 
and other methods are needed. In all cases the absolute 
values of times and coordinates, measured in the Fermi 
energy and wave length, are large. 

The rest of this paper is devoted to study of asymptotic 
properties of Calogero model in the semiclassical regime. 
We consider the Green's function first. Anyon operator 
is represented via boson operators |9j as 



$> A (x,t) ~ exp 



dx(v(x, t) + Tt\p(x, t) 



(7) 



This definition is consistent 16] with the commutation 
relations 

* A (x)* A (x') = e-^ Xs ^ x '- x ^ A (x')^ A (x). (8) 
The Green's function is given by a functional integral 



G A (x',t') 
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(9) 




FIG. 1: single particle Green's function (of right moving elec- 
trons) as a function of x for the finite band- width, v = 1, m = 
l,i= 1 
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The hydrodynamic action corresponding to the Hamilto- 
nian eq.Q is given by 



S 



dxdt 



vd x 1 d t p - 



1 
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To calculate the Green's function we use steepest decent 
method, known also as WKB. The action l|lUfl has mas- 
sive direction and low energy directions in the configura- 
tion space (a constant of shift x — > x + a,t ^ t + b) does 
not change the value of the action on any solution). The 
first may be integrated out by saddle point method, the 
latter need to be taken into account exactly. The saddle 
point equation is 

Pt + d x {pv) = -S(x - x')S(t - t') + S(x)S(t) (11) 



+ vv x + d x (^j =-Tr\S(x-x')5(t-t') + n\S(x)5(t) 

(12) 

As was suggested earlier 0, EI the elementary particle 
of fermionic theory corresponds to the quantum soliton 
of the boson theory 0] . In the limit of large A we look 
for solution (0 < t < t') in a form 

p{x,t) = l+p s {x— vt-X(t)),v(x,t)=v s (x— vt-X(t)),(13) 

where X(t) is a global variable. It corresponds to the 
motion of the soliton as a whole. This is a zero mode of 
the system. The soliton's density and velocity are given 

by m 
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Here B = ttX 2 /(v 2 — n 2 X 2 ), v = x' /t'. The soliton exists, 
provided B > 0, i.e. it's velocity exceeds the sound veloc- 
ity s = nXpo- It is very similar to the saddle point con- 
tribution for the ordinary integral of free fermions where 
saddle point contribution appeared only for supersonic 
trajectories. As velocity increases v ^> s the size of soli- 
ton B goes to zero and the particle is more localized. By 
integrating over infinitesimal region around x — 0, x = x' 
over x we see that the solution l|14|l indeed has a right 
singularity. The total number of electrons carried by the 
soliton (J p(x)dx = 1) is quantized to one. The velocity 
of soliton is not quantized. Then we exponential accu- 
racy we find that Green's function of anyon excitation in 
CS model is given by 

G A (x',t')~e iS « [ X(t) °ra e T/«'« 2 W. (15) 

JX(0)=0 

Here the value of the action on classical trajectory S c i — 
t'(v 2 /2 — sv) = '^jr — sx'. Performing the integration over 
zero mode find 



G A {x,t)~ x \- 



(16) 



where pf — npo is Fermi momentum. Equation (|16|) is 
the central result of this work. It shows that far from 
light cone the anyon's propagator oscillates as a function 
of x with the period y/t, c.f. the Green's function of free 
electrons eq. Q. 

The soliton that determines the Green's function is also 
responsible for the finite support of polarization opera- 
tor [13 H(q,t) = -i(Tp(q,t)p(-q,0)) This, as it is was 
pointed out in Ref. y], is necessary for various physical 
phenomena. Several attempts had been in the literature 
to account for it [2ll E3 | . Here we show how broaden- 
ing of the support can be understood within quantum 
hydrodynamics. 

The polarization operator is determined by neutral 
particle hole excitations. In lineal bosonization they are 
acoustic phonons. In its nonlinear version they are peri- 
odic density 0, ■ The saddle point equation on peri- 
odic waves reduces to 



A 2 
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\— + -A = C, (17) 
2 p 4 p z 



where C is an arbitrary constant. Solving this equation 
one finds 
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(18) 



Here velocity of the wave 



v = s I 1 




n B 



2B\ 
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where ub is Bose distribution function. In a momentum 
representation 



PbM = e- S|9lAr 5>M -27m). 



(20) 



From cq. !|19I20|) we find that in a long length limit 
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If the width of the picks is bigger than the distance 
between them (B 3> /) one recovers a dispersion of an- 
harmonic acooustic phonons v(q) — s — Xq/2m. This 
excitations determine the lower edge of polarization op- 
erator support. In standard bosonization only those are 
taken into account j^. 

For the waves with period longer than / 3> B there 
is only singe pick (soliton) in that length. One can as- 
sume that since other pick are far apart one can study 
its motion regardless of others. It contributes to the en- 
ergy E — £3 ^-. By recalling |^ the soliton dispersion 
(e(q) = sq + q 2 /2m, q > 0) we obtain an upper support 
of polarization operator. 
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A bosonic and a fermionic descriptions are dual 14] . 
In fermionic language excitations are characterized by 
partition k. In bosonic notations a single period neutral 
excitations are described by the parameters B and I, or 
alternatively by n and ub- The partition k = {q} corre- 
sponds to taking the topmost electrons and moving it q 
steps up and corresponds to the diluted gas of solitons. 
The conjugated partition k = \ q that describes moving 
the top q electrons by one step up corresponds to the 
anharmonic acoustic phonons. 

To summarize: the saddle point trajectory that cor- 
responds to the soliton solution of non linear hydrody- 
namic equation is dominant for anyon Green's function. 
By evaluation the value of the action on this trajectory 
we have found an explicit formula for the Green's func- 
tion of anyon. 

Due to the infinite number of periodic density waves 
a polarization operator acquires a finite support in fre- 
quency momentum plane. Of which the lower edge of 
excitations is given by dispersive acoustic phonons (such 
excitation can be accounted by linearizing an action but 
accounting phonon dispersion even for generic interac- 
tion). The upper edge is determined by a superposition 
of non interacting solitons. As the model became non 
integrable and soliton is destroyed the upper part of sup- 
ports is expected to smear. 

Note added: As this work was nearly completed such 
investigation had been performed in Ref . j2|| . It was in- 
deed found that for generic interaction the upper bound- 
ary of the support becomes smeared, while the lower 
boundary is unchanged. 

I benefited a lot the interaction with Ilya Gruzberg 
and Paul Wiegmann. I am grateful to A. Abanov, A. 
Kamenev, D. Maslov, A. Mirlin, M. Stepanov, M. Stone, 
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